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I. INTRODUCTION 



The existence of dark matter (DM) [l| provides one of the strong motivations to seek for 
physics beyond the Standard Model (SM). Weakly interacting massive particles (WIMPs) 
have been the promising candidate for dark matter, that are assumed to have weak scale 
interactions with the SM particles and a weak scale mass. However, the property and identity 
of WIMP dark matter remains unknown and is one of the big puzzles in particle physics 
and cosmology. 

Recently, it has been shown that gamma ray line spectrum exists at E^^ = 130 GeV 



using Fermi-LAT data 



groups 



0,3- 



The signature has been independently confirmed by other 



5| and officially investigated by the Fermi-LAT collaboration but with the peak 



being shifted to E^ = 135 GeV 



energetic pulsar winds [7(1 , ^ ^ 

limb signal 6|, |l3|. The Fermi-LAT collaboration [ij] 



6|. There are some possible explanations such as mono 



10 



12| including Earth 



Fermi Bubble |8|, |9| or instrumental effects 

and the H.E.S.S collaboration [1| 
have reported only the upper bound on the annihilation cross section of WIMPs, that is 



compatible with the Fermi gamma-ray line, so the dark matter interpretation of the Fermi- 
LAT line signature seems plausible. 

Dark matter is neutral and thus cannot annihilate into photons at tree-level. The gen- 
eration of photons must happen via the loops of charged particles to which dark matter is 
directly or indirectly coupled, so the annihilation cross section of dark matter into photons 
is much suppressed as compared to other tree-level annihilation channels. Thus, in order to 



realize a large branching fraction of the annihilation cross section into p 



16 



lotons, we need to 



18| or a resonance 



m. 



rely on a large coupling to new charged particles running in loops 
pole of the mediator particle between dark matter and photons 19, 

The extra f/(l) gauge symmetry is one of the simplest extensions of the SM. In this 
paper, we propose a renormalizable model of vector dark matter in which the extra U{l)x 
gauge boson couples to the SM particles through the mixing between the complex singlet 
scalar, which is responsible for f/(l)x breakdown, and the SM Higgs doublet. When there 
is a quartic coupling between the singlet scalar and a heavy charged scalar, dark matter 
can annihilate into a photon pair with sizable branching fraction, provided that the mixing 
between the singlet scalar and the SM Higgs boson is small enough. The annihilation 
cross section of dark matter into a photon pair is enhanced near the resonance pole of the 



singlet-like scalar mediator to be consistent with Fermi gamma-ray line, while the thermal 
relic density determined by the annihilation channels into WW and ZZ can be successfully 
obtained by a small Higgs mixing. We discuss various phenomenological implications of the 
model, e.g., Higgs to diphoton rate, electroweak precision measurements, collider constraints 
on the charged scalar, and the vacuum stability bound of the scalar potential. 

We note that in most of other previous works on vector WIMP, a vector boson was 



introduced as a Proca field 2l|-|23| or a scalar sector is described by the nonlinear sigma 



model (2J]. In our model, vector dark matter is realized as a gauge boson of extra U{1] 
symmetry when it is broken spontaneously by the VEV of the hidden Higgs field at the 
renormalizable level 



22 



25|. 



The paper is organized as follows. Beginning with the introduction of our model, in 
Sec. mil we calculate the annihilation cross section of the U{l)x gauge boson dark matter 
and show that the desired thermal relic density and a large annihilation cross section into 
two monochromatic photons can be realized for a consistent set of parameters. Then, we 
study the effect of the charged scalar on the Higgs diphoton decay rate in Sec. HVl and 
various experimental constraints on the charged scalar in Sec. |Vl In Sec. I VII we study the 
running of the couplings with the modified renormalization group equations (RGBs) due to 
additional interactions and consequently show how the stability of the Higgs potential is 
improved. We summarize our results in Sec. IVIII There are two appendices summarizing 
the gauge and scalar interaction vertices and RG equations, respectively. 

II. THE MODEL 

We consider a simple model of vector WIMP dark matter which couples to the SM 
particles through Higgs portal interactions. The gauge sector of the model is based on the 
SU{3)c X SU{2)l X f/(l)y X U{l)x gauge group. The extra U{l)x gauge symmetry is 
spontaneously broken by the vacuum expectation value (VEV) of a complex scalar 5*1. For 
a minimal extension of the SM ^, we introduce an SU(2) singlet charged scalar S*^, which 
carries Y = 1 but is neutral under U{l)x- We assume that all the SM particles including 
the Higgs doublet are neutral under the U{l)x- 



^ The model ean be extended with a large multiplet containing a charged scalar such as extra Higgs doublet 
or Higgs triplet, etc. 



The model has a Z2 symmetry under which S\ — )> S\ and X^ — )■ — X^, which guarantees 
the stabihty of vector dark matter X^. The Lagrangian of the model is 

£ = —F^^F^'' ^\D^Sr\^ ^\D^,S2\^ -Vi^.SuS^)^ j,,UC -L.St, (1) 

where F^y = d^X^, — d^X^, and the covariant derivatives are D^Si = {d^ — igxX^)Si and 
-0^5*2 = {dfj, — ig'Bfj_)S2, with respect to f/(l)x and f/(l)y gauge symmetry, respectively. 
After the electroweak symmetry breaking, -0^5*2 is reduced to the covariant derivative with 
respect to ?7(l)cm symmetry. The scalar potential of the SM Higgs doublet $ and complex 
scalar fields 5*1 and S2 is given by 

V{^, 81,82) = film^ + fil\8,\' + fil\82\' + h^\^\^ + ^X2\8,\^ + ^Xs\82\^ 

+A4|$n^iP + Asl^n^sp + Xg\8,\^\82\^. (2) 

Here the A4 coupling is relevant for the mixing between $ and 81, after the breaking of U{l)x 
and electroweak symmetry. The dominant annihilation of vector dark matter occurs through 
this coupling so that the correct relic density can obtained. The coupling Ag connects the 
vector dark matter to the charged particle and enhances the photon emission. The coupling 
A5 in combination of the charged scalar mass can be constrained by the branching ratio of 
Higgs boson decay into two photons at the Large Hadron Collider (LHC). 

In the last term of the Lagrangian ([T]), Lj is the SM lepton doublet with flavor index 
i = 1,2,3, C is the charge-conjugation operator and dot denotes the 8U(2) anti-symmetric 
product. This term induces the decay of heavy charged scalar 82 ^ li + Vj- The experi- 
mental constraints on the lepton flavor violating term will be discussed in the Sec |Vl This 
model can be extended to a Type-II seesaw model where the charged scalar 5*2 is embedded 



into a triplet Higgs field with Y = 2 [26|]. In this case, the lepton couphngs of the charged 
scalar would be small in order to explain neutrino masses. 

At the vacuum with a nonvanishing singlet VEV, {81) = vs/V2, the U{l)x gauge sym- 
metry is broken so the gauge boson X^ acquires mass. 

Mi = g'xvl. (3) 

We expand $ and 5*1 fields in unitary gauge, around the electroweak vacuum with ($) = 




Va/2 and (5*1) = vs/\/2, as 

/ n \ 

(4) 

(5) 

with V ~ 246 GeV. Then, two scalar modes and (ps in general mix so the mass eigenstates 
h and H are given in terms of the mixing angle a as 

/i = cosa0 — sina^s, iJ = sin a + cos a 05. (6) 

The Higgs mixing essentially depends on A4 through 



Aif2 — A2w|' 



tan2a = Y-:-^ ^— 1", (') 



and the mass eigenvalues are 

MIh = Xiv' + X2VI T ^J {\2vl - Xiv^Y + ^Xlv^vl . (8) 

For a small mixing angle a, we can regard /i ^ as being SM Higgs-like and iJ ~ 05 as 
being singlet-like, and the mass eigenvalues are Mh — Xiv'^ and Mjj ~ A2'y|. The gauge and 
scalar interactions of the Higgs fields are listed in appendix A. 

The absolute stability of the electroweak vacuum gives rise to the following conditions on 
the quartic couplings in the scalar potential J27l |. 

Ai > 0, A2 > 0, A3 > 0, 

A12 = A4 + v^V^ > 0, (9) 



Ai3 = A5 + VAiAs > 0, 



A23 = As + V A2A3 > 0, 
and 

-^123 = -y ^1^2-^3 + ^6V '^l + ^5V ^^2 + '^4V -^3 + V -^12-^23-^13 > 0. (10) 

Throughout the work, we also impose perturbativity conditions until the Planck scale as 

|A,|<47r. (11) 
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FIG. 1: Feynman diagram for annihilations of vector dark matter into WW^ ZZ at tree level. 

III. RELIC DENSITY AND FERMI GAMMA-RAY LINES 

In this section, we compute the DM annihilation cross sections and discuss the constraints 
of the relic density and the Fermi gamma-ray line in the model. 

A. Thermal relic abundance 



The thermal relic abundance of the vector WIMP dark matter X is estimated by inte- 
grating the following Boltzmann equation for the dark matter number density nx in the 
early Universe 28|, |29| , 

dnx 



dt 



3Hnx = -{av)(nj,-iniy^ 



(12) 



where H, {av), and Uj^ denote the Hubble parameter, the thermal averaged annihilation 
cross section times relative velocity, and the dark matter number density at thermal equi- 
librium, respectively. X dominantly annihilates into W and Z boson pairs through the 
s-channel exchange of the Higgs bosons h and H, as shown in Fig. [1] The magnitude of 
those annihilation cross sections is proportional to (sin a cos a)^ and thus scaled by the mix- 
ing angle between Higgs bosons. A larger (smaller) annihilation cross section is realized for 
a larger (smaller) sin a. The velocity times DM annihilation cross sections into W, Z boson 
pairs before thermal average are given by 

(13) 
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FIG. 2: The plot of relic density of the dark matter Qh^ vs. its mass. Here we used Mg± = 140 GeV, 
vs = 1 TeV, sin a = 0.00047 and Ag = 0.2. We can ignore the A5 dependence for the parameter 
region in our interest. The horizontal line is the relic density of cold dark matter revealed by the 
Planck result, Ocdm/i^ = 0.1199 ± 0.0027 [SO]. 
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QxxhQhzz 



Ml + IMhTh 
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(14) 



where s is the total energy at the center of mass frame and each coupUngs are given by 



9xxh = -'^gx^'s sin a, 
1 o 



9hWW 

Qhzz 



-g2Vcosa, 



9xxH = '^g.x^'s cos a, 



9HWW = 7Tfi'2^'sma, 



(15) 
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2 cos 9ly 



V cos a, 



9hzz 



9l 



2 cos O^Y 



f sma. 



The heavy Higgs boson mainly decays into XX, and also into W, Z and h pairs suppressed 
by the mixing. The decay width is given by 
,2 
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FIG. 3: Feynman diagrams for the annihilation of vector dark matter into a photon pair at one-loop. 



We are interested in the case where X has a large annihilation cross section into two 
photons so that the gamma ray line at 135 GeV observed by Fermi LAT can be explained. 
Taking the annihilation into a.W or Z boson pair through h and if -exchange to be strongly 
suppressed due to a tiny mixing angle between SM Higgs boson and singlet scalars, namely, 
I sin a I ^ 1, we obtain the desired thermal relic density and the necessary cross section 
into a photon pair for Fermi gamma-ray line near the resonance pole of the singlet-like 



scalar 



19 



20| . Therefore, near the resonance, Mx — Mh/2, for the WIMP mass^ around 



135 GeV, we need to take the singlet- scalar mass M^ to be about 270 GeV. 

In Fig. m we show that the right relic density of WIMP can be obtained by performing 
the thermal average of the annihilation cross section with the procedure in Ref. 29|. We 
note that there is no bound on vector dark matter from direct detection in our model, due 
to a tiny Higgs mixing. 



B. Monochromatic photons from DM annihilation 



Now we discuss the aspect of the indirect detection of vector WIMP. As mentioned in 
the previous subsection, X annihilates mainly into W and Z boson pairs. In addition to 
those modes, X can also annihilate into two monochromatic photons with the effective 
interaction, which is induced by loops with the charged scalar S*^ as shown in Fig. [31 Hence, 
this diphoton mode can provide a source for the gamma-ray line observed by Fermi LAT. 
When the singlet-like scalar H has small couplings to W and Z bosons due to a tiny sin a 
but it has a sizable coupling to the charged scalar 5*2, the diphoton mode takes a larger 
branching fraction of the annihilation cross section than usually expected. 



^ Of course, there is another resonance region Mx — AIii/2 ~ 63 GeV where the main annihilation mode 
is XX -^ bb. 




FIG. 4: Thermal averaged annihilation cross section (Dashed line), given by the sum of 77, WW, ZZ 
final states at the freeze-out temperature of dark matter, Tj = Afx/20. Solid lines are at zero 
temperature (Red: 77, Blue: WW, Green: ZZ final states). Here we used the same parameters in 
Fig. [2j The horizontal black line is the value of the DM annihilation cross section into two photons, 
(^'^)xx-5'77 ~ ^-^ ^ 10"^*^ GeV~ , needed for the Fermi-LAT gamma ray line. 

The velocity times DM annihilation cross section into a photon pair before thermal av- 
erage is expressed by 



{av)^^{s) 



Otr, 



QGtt^s 
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s-Mi__+iMHTH, 



Ml 



'Ml 



1 - 
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2Ml V2 
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(17) 



with 



/(r) = [sin-^(l/v/7)]^ 



(18) 



Keeping only the //-exchange contribution in Eq. ( TT7|) for a small Higgs mixing and using 
s ^ 4:Ml, we get an approximate form for the thermal averaged annihilation cross section 
into a photon pair at present as 

2 



{(Tv) 



at. 
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967r3M|- 



Xq cos^ aMl 



AMI -Mh + iMnTH 



Mi 
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Ml 



■f 



Ml 



(19) 



In Fig. m we show, as a function of the DM mass, the annihilation cross sections into 
WW, ZZ, '-ff at zero temperature, and the thermal averaged total annihilation cross section 
used to estimate Qh"^, respectively. To explain the gamma ray line spectrum of the Fermi- 
LAT 2|, |3|, for the Einasto dark matter profile, we require that 

{av)^^^^^ = (1.27±0.32+°i^) x 10-2^cm=^s-\ 
^ 1.1 X 10"^°GeV~l 

In order to obtain the observed Fermi gamma-ray line together with the correct relic 
density, the small mixing angle, | sina| <^ 1, is necessary. Furthermore, the annihilation into 
WW and ZZ modes have to be suppressed enough not to generate too many continuum 



photons 



31 



1|-|33| . In our case, the DM annihilation cross sections into WW and ZZ are 
suppressed at present for the parameters which explain both the Fermi gamma-ray line and 
the relic density. For instance, for Mx = 134.74 GeV, we obtain (cv)^^ = 1.09 x 10~^''cm^/s 
and Qh^ = 0.118 while {<yv)zz/ {o'v)^^ = 0.43 and {av)v/Y/ / {av)^^ = 0.28. We note that the 
total DM annihilation cross section at present is smaller than thermal cross section, because 
the temperature effect shifts the peak of the resonance at freezeout towards a smaller DM 



mass as compared to the case with zero temperature |34l ]. 

The same diagrams in Fig. [3] applies to the annihilation of vector dark matter into Z'j 
and loop induced ZZ final states. The annihilation XX — )■ Z'-f emits additional gamma 
line at E^ = 114 GeV and the resulting flux is suppressed by 0.21 as compared to that of 
130 GeV gamma ray line [l7|, while the loop induced annihilation into ZZ is negligible. 

IV. HIGGS TO DIPHOTON RATE 

The charged scalar 5*2, introduced to explain the Fermi gamma-ray line, can give a positive 
or negative sizable contribution to the Higgs to diphoton rate due to the quartic coupling 
A5 to the SM Higgs field. In this section, we discuss the constraint on the modified decay 
rate of Higgs boson h to diphotons from the recent measurements at the LHC. 

We define the ratio of the Higgs production cross section times the branching fraction, 
/i^^ = , ^^^ ''Y — • The reported signal strengths for the Higgs to diphoton rate from ATLAS 



and Multi-Variate- Analysis (MVA) of CMS data are the following 35 1 



ATLAS _ 1 «c;+0.34 CMS _ q 70+O.28 /qi N 

P77 —-'-•'-'^-0.30' P77 — U.^O_o.26- \^^) 
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FIG. 5: Contours of Higgs to diphoton rate on A5 — M^± plane. The green and yellow regions 
correspond to the combined LHC bounds on R^^ at 68% and 90% C. L., respectively. 



Following a similar method as in Refs. [36|, |37[ and assuming that the combined data is 
Gaussian, we have derived the combined value of the Higgs to diphoton rate as 



/i™'^^' = 1.18 ±0.20. 



The SM-like Higgs boson decay width F(/i — )■ 77) is given by 



38 



V{h ^ 77) 



128V27r3 



gwAi{Tw) + gtQtNcAii2{Tt) + ghAo{Tg±] 



with loop functions 



-x2[2x-2 + 3x-^ + 3(2x-^ - l)/(x~^)], 



Ai(x) 
^1/2(3;) 
Ao{x) = -x^[x~^ - f{x~^)] 



fix) 



arcsin^ y/x for x < 1 



-i(ln((l + VI - .^"^)/(l + VI - x-i)) - i-nf fora;>l, 



(22) 

(23) 

(24) 
(25) 
(26) 

(27) 
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and Ti = AM^/Ml. Qt = ^, Nc = 3 ioi top quark, gw and g^ are the Higgs trilinear couplings 
to W gauge boson and top quark normalized to the ones of the SM, and in our case those 
are almost 1. The Higgs coupling to the charged scalar boson is 



M, 



w 



9H = -j^vX,. (28) 



gM 



.± 



By taking the ratio to the SM value, we obtain 






'' - Tih ^ 77) 



SM 



^ 9hAo{Ts±] 



gwAiirw) + gtQ^tN.Ay^irt) 



(29) 



which is a function of A5 and Mg±. When the production cross section of h is SM Higgs- 
like, the above ratio is approximated to the signal strength measured at the LHC, that is, 

In Fig. \5\ we depict the parameter space (A5, M^±), showing the contours of the h decay 
to diphoton rate with other Higgs couplings being assumed the same as in the SM. We find 
that, from the combined value of ATLAS and CMS (MVA) diphoton signal strengths at 90% 
C.L., the extra quartic coupling is constrained to —2.5 ^ A5 < 0.7 for Mg± = 140 GeV. The 
heavier the charged scalar, the larger values of the extra quartic coupling A5 are allowed. 

V. CONSTRAINTS ON CHARGED SCALAR 

We have introduced the lepton couplings to the charged scalar 5*2 by gauge invariant 
terms, fijL^ ■ LjS2, so 5*2 could be unstable. In this section, we discuss the phenomenology 
of the charged scalar from the indirect limits and the collider search for charged particles at 
the LHC. 

The lepton couplings to the charged scalar are similar to lepton number (R-parity) violat- 
ing terms in the minimal supersymmetric standard model, AW = XijkLi ■ LjE^, so the same 
bounds from precision measurements are applied to them. The bounds from the charged 
current universality are I/12I < 0.04(?72^±/(100GeV)), and the constraints from Rr = r(r — >■ 
ei^z/)/r(r — )■ fiuu) or i?,-^ = T{t — )■ jiui))/V{ii — )■ euu) are \fij\ < 0.05(mg±/(100 GeV)), and 
Vii deep inelastic scattering gives the bound, I/12I < 0.02(mg±/(100 GeV)) |39(]. Other lepton 



M 



flavor violation Br{^ — )■ 67) also gives a similar bounds 40|]. The lepton Yukawa couplings 



12 



gives a negative contribution^ to the muon anomalous magnetic moment as 4ll. |42| 



?Ti 1 

Then, using the bounds from precision measurements, we get a very small contribution, 
|Aa^| < 3.45 X 10^^^. We note that as far as the electroweak precision bounds on fij 
are satisfied for Mg± ^130 GeV, the continuum photons coming from the three-body DM 
annihilation into 5*2/// can be suppressed enough. 

New particles with electroweak charges have been searched for at colliders. The strin- 
gent bounds on the charged scalar come from the direct slepton pair production at the 
LHC 



43 



44| . where a left-handed slepton can decay into lepton and neutralino. The 
opposite-sign dilepton search with the same-flavor channel at CMS excludes slepton masses 
between 110 GeV and 275 GeV for massless neutralino 4^. But, in our case, the charged 



scalar can decay into all the charged leptons: S2 -^ eUf^^r, fJ''^e,T, and Ti?e,fi- Since the lepton 
coupling matrix, fy, is anti-symmetric, at least two different flavors always appear in the 
decay product of the charged scalar. Therefore, the CMS mass limit with the same-flavor 
channel scales down or does not apply, depending on the branching fraction of the same- 
flavor decay mode. Currently, the most stringent constraint on the charged scalar mass 



comes from the LEP exclusion limit up to 95 GeV 45| 



VI. VACUUM STABILITY 



The discovered scalar boson with 126 GeV mass has been shown to have very similar 



properties to the SM Higgs boson with more precision 46|. Although we need more data 
to confirm the properties of the Higgs boson, we assume that the discovered scalar boson is 
SM Higgs-like. In this case, a small Higgs quartic coupling leads to a problem of vacuum 
instability below the Planck scale J47| , requiring new physics beyond the SM. In this section, 
we discuss the effect of the additional quartic couplings between the Higgs boson and extra 
scalars in the model, taking account of dark matter constraints from Fermi gamma-ray line, 
Higgs boson data and other collider bounds, discussed in the previous sections. 



^ There was an error in the previous works on the LLS2 coupUng [40[ which showed a positive contributii 
to the muon anomalous magnetic moment. 
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FIG. 6: Running couplings as a function of t = ln(^/Aff). We have chosen Ai = 0.26, A2 = 
0.07, A3 = 0.15, A4 = -0.0001, A5 = 0.28, Ae = 0.2 and gx = 0.134 at Mt, which leads to Mx = 
134 GeV, Mh = 270 GeV, Mh = 126 GeV and sin a = 0.0005 with vs = ITeV. The lepton 
couplings to the charged scalar are ignored in the RG analysis. 



e s. 



In our model, it is possible to have a sizab 
presence of the mixing with a singlet scalar 48N50| 



lift in the Higgs quartic coupling in the 
as follows. 



A, 



cff 



\i-6\ 



with 



2\2 



6X 



{Ml - Ml) 



sm a cos a 



(31) 



(32) 



v'^^Mfj cos2 a + Ml sin^ a) ' 
In the decoupling limit of a heavy singlet-like scalar with Mh ^ Mh-, the tree- level shift 
is approximated to 5\ ~ M|f sin^ a/(f^) ~ AI/A2 |48|, |49[. In this paper, however, we take 
the singlet-like scalar mass to be close to the resonance pole, Mh ~ 2Mx ~ 270 GeV. 
Furthermore, since | sina| <^ 1 for the correct relic density at the resonance, the tree- level 
shift in our case becomes 5X ~ 0.7 sin^ a, which is extremely small. 

Now we consider the RG effect on the Higgs quartic coupling. As shown in Appendix [Bl 
there are positive contributions to the beta function of the Higgs quartic couphng, A4 and 
A5, in the RG equations, so the vacuum stability can be improved as compared to the SM. 
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FIG. 7: Vacuum stability conditions: A12 > 0, A13 > 0, A23 > and A123 > 0. Tlie same parameters 
are chosen as in Fig. [6l 

But, from the results of the previous sections, the quartic couphng A4 between the Higgs 
doublet and the singlet scalar must be small because of the relic density condition, hence its 
contribution to the running of the Higgs quartic coupling is negligible. On the other hand, 
a sizable quartic coupling A5 between the Higgs doublet and the charged scalar is allowed, 
being consistent with the observed Higgs boson decay rate to diphoton. 

If A5 is positive, it can help increase the vacuum instability scale, without violating the 
new vacuum stability conditions of extra scalars. We note, however, that if A5 is negative 
and large as suggested by Higgs data, it could increase the Higgs quartic coupling by the 
RG further, but perturbativity bound and extra vacuum stability conditions strongly restrict 
this possibility. In Fig. El we show the running couplings until the Planck scale for the low- 
energy couplings including a positive A5, that are consistent with Fermi gamma-ray line, 
relic density, Higgs diphoton data and indirect and collider bounds. In Fig. [TJ the vacuum 
stability conditions are shown to be satisfied until the Planck scale, for the same parameters 
as in Fig. [61 
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VII. CONCLUSION 

We have proposed a renormalizable model of vector dark matter where the extra U(l)x 
gauge boson is a dark matter candidate and it interacts with the SM particles through the 
Higgs portal term, namely, the mixing between the U{l)x breaking singlet scalar and the 
SM Higgs doublet. If the Higgs mixing is small enough, the DM annihilations into W and Z 
boson pairs at the resonance pole of the singlet-like scalar can reproduce the correct thermal 
relic density without overproducing continuum photons. In the presence of a quartic coupling 
between the singlet scalar and the charged scalar 5*2, vector dark matter also annihilates 
into a photon pair with a sizable branching fraction at the same singlet resonance. 

As long as the couplings of the charged scalar to the SM leptons are small enough, i.e. 
fij < O(10~^), we showed that all the electroweak precision constraints concerning leptons 
are satisfied. Even though it would be very difficult to find the singlet scalar at colliders 
due to a tiny mixing with the Higgs, the charged scalar would be accessible at the LHC or 
linear colliders, due to a distinct signature that two opposite-sign leptons of different flavors 
are equally produced from the charged scalar decay. We have also shown that the vacuum 
stability bounds are satisfled until the Planck scale, due to a sizable Higgs coupling to the 
charged scalar, which is allowed by the current Higgs diphoton data. 
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Appendix A: Scalar interaction vertices 

Gauge interactions: 

^int = gxX^vs{H COS a — hsina) 



+ j^{2clrW-W^ + Z^)v{hcosa + Hsina). (Al) 



Scalar interactions: 



—Cint =-{vXi cos^ a + f A4 sin^ a cos a — vs\2 sin'^ a — vs\i sin a cos^ Q.)h^ 

H — (t;sina;(3Ai + A4) + 3f sin3a(Ai — A4) + f5Cosa(3A2 + A4) + 3f5Cos3a;(A4 — \2))h?'H 
8 

+ -(i;cosa(3Ai + A4) — 3t'cos3tt(Ai — A4) — fs'Sina(3A2 + A4) + 3f5sin3a(A4 — \2))hH^ 
8 

H — (f Al sin'^ a + wA4sinacos^ a + vs\2 cos^ a + f5A4Cosasin^ a)H^ 

+ {v\^ cos a — vs\q sin a)hS^S^ + {vX^ sin a + t'^Ag cos a)HS^S^ 

H — (cos^ aAi + sin^ aA2 + 2 cos^ a sin^ aA4)/i^ H — (sin^ aAi + cos^ aA2 + 2 cos^ a sin^ aA4)if ^ 
8 8 

H — (cos"^ asina(Ai — A4) — cosasin'^a(A2 — \4))h^H 

+ -(cosasin^a(Ai — A4) — cos'^asina(A2 — \4))hH^ 

-\ — ((cos^a — 4cos^ asin^ a + sin^ a)A4 + 3cos^ asin^a(Ai + \2))h?H'^ 

+ -(cos^ aAs + sin^ a\Q)h?' S^ S~ + -(sin^ a\^ + cos^ a\Q)H'^S^ S^ 

+ sin a cos a (A5 — \)hH S^ S~ + -AajS'^^S'^p. 

(A2) 

Appendix B: Renormalization group equations 

The running of the couphng pj is governed by the RG equation, -^ = /3p. with /3p. being 
the corresponding beta function and t = ln(/i/mi). The beta functions for scalar quartic 



17 



couplings with a = IGtt^ are 

«:/3,, = 12A? + {12y^ - 9/ - 3g")\^ - Uyf + ^-g' + \g'' + ^^V + 2A^ + 2A^, (Bl) 
«:/3a, = 18A^ + \\\ + 2A2 - 12(7iA2 + 12(7^, (B2) 

«/3a3 = 18A^ + 4A^ + 2A^ + (4Tr(/t/) - 12(?'2) A3 + 12^7'^ - 4Tr(//t//t), (B3) 

«:/3a, = (12Ai + 6A2)A4 + 2A5A6 + (^vl - ^g'^ - ^g' - 6(7^) A4 + g'^gl, (B4) 

kP,, = (12Ai + 6A3)A5 + 2A4A6 + {6y^ + 2Tr(/t/) - ^g'^ - ^g^^, + 9g'' + 9(?VtP5) 

«:/3ac = 6(A2 + A3) As + 4A4A5 + (2Tr(/t/) - 6g" - 6gj,^ Ag + Ug'^gj,. (B6) 

The beta functions for the top Yukawa couphng and the lepton Yukawa couphngs to the 
charged scalar are 

-/?.. = yt{ly^-89l-l9'-^^9"), (B7) 

«:% = 4(//t/)., + /., (4Tr(/t/) - ^/ - ^'^ . (B8) 

Here, we have ignored the charged lepton Yukawa couplings to the SM Higgs field. When 
a single lepton coupling to the charged scalar, e.g. / = I/12I, is dominant, the RG equation 
for that becomes 

^Pf = f{l2f-lg'-lg"). (B9) 

The beta functions for the gauge couplings are 

43 19 1 

^/3g' = Y^'^ K(3g = -Y^^ f^(3gs = -7gl K(3g^ = -g^x. (BIO) 
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